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ABSTRACT 

Context. V1094 Tau is bright eclipsing binary star with an orbital period close to 9 days that contains two stars similar to the Sun. 
Aims. Our aim is to test models of Sun-like stars using precise and accurate mass and radius measurements for both stars in V1094 Tau. 
Methods. We present new spectroscopy of V1094 Tau, which we use to estimate the effective temperatures of both stars and to refine 
their spectroscopic orbits. We also present new, high-quality photometry covering both eclipses of V1094 Tau in the Strdmgren uvby 
system and in the Johnson V-band. 

Results. The masses, radii, and effective temperatures of the stars in V1094 Tau are found to be = 1.0965 ± 0.0040 Afg, = 
1.4109 ± 0.0058 Ro, Teff,A = 5850 ± 100 K, Mg = 1.0120 ± 0.0028 Mq, Rb = 1.1063 ± 0.0066 Rq, and Teff,B = 5700 ± 100 K. An 
analysis of the times of mid-eclipse and the radial velocity data reveals apsidal motion with a period of 14500 ± 3700 years. 
Conclusions. The observed masses, radii, and effective temperatures are consistent with stellar models for an age 6 Gyr if the stars 
are assumed to have a metallicity similar to the Sun. This estimate is in reasonable agreement with our estimate of the metallicity 
derived using Stromgren photometry and treating the binary as a single star ([Fe/H]= -0.09 ±0.11). The rotation velocities of the stars 
suggest that VI094 Tau is close to the limit at which tidal interactions between the stars force them to rotate pseudo-synchronously 
with the orbital motion. 
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1. Introduction 


VI094 Tau (HD 284195) is a ninth-magnit ude star that w as dis¬ 
covered to be an eclipsing binary star by iKaiserl d 19941) using 
photographic photometry. The period quoted in that paper is in¬ 
correct and was later found to be the time between the secondary 
and pri mary minima. The corr ect period {P ^ 8.988 d) was first 
given bv iKaiser & Frevid 19981) . who also established that the or¬ 
bit is eccentric, with the secondary minimum falling near phase 
0.65. Observations of the times of mid-ecli pse have continue d 
since its discovery and have been analysed bv lWolf et al.l d2010h . 
who claim to have detected apsidal motion with a period of 
13600 years and variations in these times of mid-eclipse with 
a period of about eight years and an amplitude of 200 s that they 
suggest may be due to a third body in the system. 

The spectrosc opic orbits of both stars have been measured to 
good accuracy bv iGriffin & BoffinI d2003h . No good light curves 


* The light curves shown in Fig. 1 are available in electronic form at 
the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A± A/ 

** Deceased 2011 June 5 


of the star were available at that time, but they were able to 
establish that the rotation of the stars is likely to be pseudo- 
synchronised with the orbital motion, that the masses are close 
to 1.1 Mq and I.QMq, and that the spectral types of the stars are 
approximately GO and G3. 

The techniques for measuring the masses and radii of the 
stars in eclipsing binaries, such as VI094 Tau to good accuracy 
(ss 1%) independently of any stellar structure models, are now 
well established. To obtain masses and radii to this precision re¬ 
quires light curves of good quality that cover both eclipses and 
spectroscopic orbits with good phase coverage based on high- 
resolution spectroscopy. With data of this quality, it is possi¬ 
ble to critically test stellar structure models, particularly if these 
mass and radius estimates can be complemented with reliable 
estimates f or the effective tem peratures of the stars and their 
metallicity (iTorres et al.ll2010h . Compilations of such data have 
also been used to establish empirical relationships between mass, 
radius, effective temperature, etc. These empirical relationships 
can be used, for example, to estimate the mass of a star based 
on its observed density, effective temperature, and metallicity. 
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These quantities a re directly measur able for stars in transiting 
exoplanet systems (lEnoch et al.ll2010|) . 

In this work we present the first high-quality light curves and 
the first high-resolution spectroscopy for VI094 Tau. We use 
these data to measure the masses and radii of the stars in this 
binary to better than 1%. We also make useful estimates of the 
effective temperatures of the stars and their metallicity based on 
Stromgren photometry. 

2. Observations 

2.1. Spectroscopy 

VI094 Tau was observed with a Cassegrain-mounted echelle 
spectrograph attached to the 1.5-m Wyeth reflector at the Oak 
Ridge Observatory (Harvard, Massachusetts, USA). A total of 
59 useful exposures were obtained between 1995 December and 
1997 April at a resolving power of R 35,000. A single echelle 
order was recorded with an intensified photon-counting Reti- 
con detector, giving 45 A of coverage centred at 5187A, and in¬ 
cluding the lines of the Mgib triplet. The signal-to-noise ratios 
range from 13 to 24 per resolution element of 8.5kms '. The 
wavelength calibration was established by means of exposures 
of a Th-Ar lamp before and after each science exposure, and all 
reductions were carried out with st andard proced ures as imple¬ 
mented in a dedicated pipeline (see lLathamlll992h . 

Spectroscopic observations were also obtained in 2002 Octo¬ 
ber using the 2.5-m Isaac Newton Telescope (INT) on La Palma. 
The 500 mm camera of the Intermediate Dispersion Spectro¬ 
graph (IDS) was equipped with a holographic 2400 lines mm ' 
grating. An LEV 4kx2k CCD was used and exposure times 
were 300 s. Erom measurements of the full width at half max¬ 
imum (EWHM) of arc lines taken for wavelength calibration, 
we estimated that the resolution is 0.2 A. A total of 64 spectra 
were taken covering the interval 4230^500 A, with estimated 
signal-to-noise ratios of approximately 50 per pixel. 

The re duction of all spectra was undertaken using optimal 
extraction (lHorndT986h as imp lemented in the software tools 
PAMELA and MoujJ (!Marshl[l989h . 

2.2. Photometry 

The differential uvby light curves of V1094 Tau were observed 
at the Stromgren Automatic Telescope (SAT) at ESO, La Silla 
and its six-channel uvbyjS photometer on 72 nights between 
October 2000 and January 2008. They contain 670 points per 
band. HD 26736 and HD 26874 were used as comparison stars 
throughout. In addition HD 27989AB was used as comparison 
star until JD2451889, but was found to be variable and was 
therefore replaced b y HD 24702. HD 27989 AB is now known 
to be a BY Dra star Istrassmeier et al.l|2000t) . HD 24702 is con¬ 
stant within the observational accuracy, whereas HD 26736 and 
HD 26874 scatter slightly more than expected; see Table [1] Eor 
HD 26874 the effect seems to be random, whereas HD 26736 
brightened by about 0.02 mag during one observing season. The 
light curves are calculated relative to HD 26874. Observations 
of HD 26736 (except for the bright period) and HD 24702 (when 
observed) were also used, shifting them first to the light level of 
HD 26874. 

We also obtained light curves of VI094 Tau with the Web- 
Scope instrument at the NE/ Observatory (NEO) located near 

’ PAMELA and MOLLY Were written by Prof. Tom Marsh and are available 
at www.warwick.ac.uk/go/trmarsh 


Silver City, New Mexico and the Undergraduate Research Stud¬ 
ies in Astronomy (URSA) WebScope on the roof of the Kimpel 
Hall on the University of Arkansas campus at Eayetteville. The 
NEO instrument consists of a 24-inch Cassegrain reflector with 
a field-widening c orrecting lens near the focus housed in a roll¬ 
off roof structure (iGrauer et al.ll200^ . At the focus is a camera 
based on the Kodak KAE-4301E charge coupled device (CCD) 
with a field of view of about 27x27 arcmin. VI094 Tau was ob¬ 
served at the NEO on 116 nights between 2 Jan 2006 and 24 Mar 
2012, producing a total of 5714 observations from 30 second ex¬ 
posures with a Bessel V filter. 

The URSA WebScope uses a 10-inch Meade LX200 
Schmidt-Cassegrain telescope scope with a Santa Barbara Imag¬ 
ing Group (SBIG) STS CCD camera, housed in a Technical In¬ 
novations RoboDome enclosure. VI094 Tau was observed on 
91 nights between 6 Mar 2001 and 25 Mar 2012, producing a 
total of 8085 observations from 30 second exposures, also with 
a Bessel V filter. 

The images were automatically measured by using the appli¬ 
cations Multi-Measure and Measure written by author Lacy. 
The software was used to locate the stars of interest in the cali¬ 
brated images and to perform background subtraction and aper¬ 
ture photometry in a region 22 arcsec square around each star 
in the NEO images and 30 arcsec square for the URSA images. 
The comparison star used was HD 284196 and the check star 
was HD 284197. Differences in the atmospheric extinction were 
corrected for all three stars. To form the differential magnitudes 
for analysis we used the flux of both comparison and check star 
for the NEO images, whereas in the URSA images only the com¬ 
parison star HD 284196 was used. The differences between the 
comparison and check star magnitudes averaged 0.008 mag for 
the URSA images on 88 nights and 0.006 mag for NEO images 
on 118 nights. These magnitude differences are similar to the 
residuals of the light curve model fits derived from the analysis 
below. The complete light curves obtain using the NEO, URSA 
and SAT telescopes are shown in Eig. [T] 

3. Analysis 

3.1. Eclipse ephemerides 

We used the jkteboi0 light curve model (ISouthworthll2013() to 
measure times of mid-eclipse from our SAT, NEO, and URSA 
photometry. We identified nights of data where our observa¬ 
tions covered the majority of the either the primary or secondary 
eclipse. All the parameters of the light curve model except the 
time of mid-eclipse and the zero-point of the magnitude scale 
were fixed at values similar to those given in Section [3. 3. 3I for 
our adopted light curve solution. We then used a least-squares fit 
to the data from each night individually to determine the times of 
mid-eclipse given. The standard errors on these times were es¬ 
timated using the cyclic residual permutation method (“prayer- 
bead” method). 

We used a least-squares fit of a linear ephemeris to these new 
times of mid-eclipse and previously published values to obtain 
the optimum values for the reference times of primary and sec¬ 
ondary eclipse and two independent estimates of the orbital pe¬ 
riod, one from the primary eclipses and one from the secondary 
eclipses. If we use the standard errors quoted on the published 
times of minimum and from the prayer-bead method for our new 
times of minimum we find that the reduced chi-squared value for 
the fit is much greater than 1. Our assumption is that this is due 

^ WWW.astro.keele.ac.uk/~jkt/codes/jktebop.html 
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Table 1. Photometric data for V1094Tau and the comparison stars. For V1094Tau, the uubyfi information is the mean value outside eclipses. N 
is the total n umber of obse rvations used to form the mean values, and tr is the rms error (per observation) in mmag. References are: M15 = This 
paper, 094 = l01senl dl994l) . 


Object 

Sp. Type 

Ref. 

V 

cr 

b-y 

cr 

m\ 

cr 

Cl 

cr 

N{uvby) 


cr 


VI094 Tau 

GO 

M15 

9.020 

9 

0.415 

4 

0.199 

9 

0.330 

11 

206 

2.596 

5 

12 

HD 26736 

G5 

M15 

8.050 

8 

0.408 

5 

0.235 

8 

0.344 

5 

113 

2.602 

5 

13 



094 

8.034 

5 

0.414 

3 

0.223 

4 

0.332 

6 

1 




HD 26874 

G4 V 

M15 

7.824 

8 

0.432 

4 

0.270 

7 

0.315 

5 

171 

2.596 

5 

11 



094 



0.438 

3 

0.269 

4 

0.320 

6 

1 




HD 27989AB 

G3V-HG6V 

M15 

7.517 

13 

0.424 

4 

0.253 

6 

0.329 

4 

67 

2.596 

6 

16 



094 

7.532 

5 

0.424 

3 

0.246 

4 

0.326 

6 

1 




HD 24702 

GO 

M15 

7.844 

5 

0.422 

5 

0.234 

9 

0.352 

5 

48 






094 



0.435 

3 

0.225 

4 

0.354 

6 

1 






Phase from primary minimum 


Fig. 1. Light curves of V1094 Tau obtained with the SAT (Strbmgren u, v, b, and y) and at the NFO and URSA observatories (V-band). For clarity, 
differential magnitudes are plotted offset in multiples of 0.5 magnitudes. 


to systematic noise sources that are both instrumental and astro- 
physical in origin, and that all the times of minimum are equally 
affected by these noise sources. In order to achieve a reduced 
chi-squared value xf ~ ^ for the least-squares ht to the times of 
primary eclipse we have added 59 s in quadrature to the standard 
error estimates quoted on published values and also to the esti¬ 
mated standard errors for our new times of minimum. For the 
times of secondary eclipse we include an additional 52 s in the 
standard error estimates. The data used and the adopted standard 
errors are given in Table |2] The resulting linear ephemerides for 


the heliocentric Julian Dates (HJD) of primary and secondary 
eclipse, respectively, are 

HJDCTpri) = 2454699.33513(19)-H 8.9885445(12)£ 
and 

HJDCTsec) = 2454498.46151(23)-H 8.9885515(13)£, 

where E is the cycle number and the hgures in parentheses are 
the standard errors in the two final digits for each free parameter. 
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Table 2. Times of mid-eclipse for VI094 Tau. Times are given as UTC 
heliocentric Julian date (HID) and are labelled as primary (“P”) or sec¬ 
ondary (“S”) eclipses. Note that the adopted standard errors used here 
are not those quoted in the original sources. 


HJD-2400000 


Type 

Source 

49653.632 

+ 0.003 

S 

Kaiser et al. 119951 

49656.762 

+ 0.003 

P 

Kaiser & Frev 119981 

49680.597 

+ 0.003 

S 

Kaiser et al. (19951 

49683.727 

+ 0.003 

P 

Kaiser & Frev 119981 

49701.7061 

+ 0.0014 

P 

Kaiser & Frev 119981 

49707.5649 

+ 0.0019 

S 

Kaiser et al. (19951 

49710.6950 

+ 0.0021 

P 

Kaiser & Frev 119981 

49755.6355 

+ 0.0021 

P 

Kaiser & Frev 119981 

50456.7393 

+ 0.0021 

P 

Kaiser & Frev 119981 

50474.7156 

+ 0.0021 

P 

Kaiser & Frev 119981 

52278.2875 

+ 0.0011 

s 

Wolfet al. 120101 

52362.3148 

+ 0.0011 

p 

Wolfet al. 120101 

52601.87780 

+ 0.00140 

s 

V, URSA 

52628.84325 

+ 0.00062 

s 

V, URSA 

52637.83158 

+ 0.00065 

s 

V, URSA 

52898.4989 

+ 0.0011 

s 

Wolfet al. 120101 

52997.3708 

+ 0.0028 

s 

Hubscher 120051 

53045.4406 

+ 0.0035 

p 

Hubscher et a". 120051 

54438.66519 

+ 0.00079 

p 

V, URSA 

54447.65487 

+ 0.00098 

p 

V,NFO 

54447.65613 

+ 0.00062 

p 

u, SAT 

54447.65604 

+ 0.00073 

p 

V, SAT 

54447.65602 

+ 0.00093 

p 

b, SAT 

54447.65617 

+ 0.00130 

p 

y, SAT 

55157.75013 

+ 0.00090 

p 

V,NFO 

55175.72858 

+ 0.00062 

p 

V,NFO 

55181.59081 

+ 0.00071 

s 

V, URSA 

55202.69394 

+ 0.00086 

p 

V,NFO 

55208.55713 

+ 0.00069 

s 

V, URSA 

55247.63640 

+ 0.00086 

p 

V,NFO 

55274.60204 

+ 0.00058 

p 

V, URSA 

55813.91459 

+ 0.00101 

p 

V, URSA 

55831.89207 

+ 0.00103 

p 

V, URSA 

55831.89275 

+ 0.00073 

p 

V,NFO 

55849.8700 

+ 0.0015 

p 

Diethelm (2012al 

55855.73235 

+ 0.00079 

s 

V, URSA 

55882.69860 

+ 0.00057 

s 

V,NFO 

55894.81080 

+ 0.00070 

p 

V, URSA 

55936.62917 

+ 0.00085 

s 

V,NFO 

55945.61739 

+ 0.00078 

s 

V,NFO 

55945.62050 

+ 0.00080 

s 

Diethelm 1201261 


The difference in the orbital period obtained from the times 
of primary and secondary eclipse is significant (4-(r) and is char¬ 
acteristic of binary stars in which the orientiation of the eccentric 
orbit is changing (apsidal motion). 

We also analysed the photometry of VI094 Tau described by 
IWraight et al.l (1201 ll) but found that these data were too noisy to 
add any useful constraints on the ephemeris or the radii of the 
stars so we do not discuss them further here. We see no evidence 
for the periodic variation in eclipse times claimed by I Wolf et al.l 

(I 2 OIOI 1 . 


Table 3. Radial velocities of V1094 Tau from CfA spectroscopy. 


HJD-2400000 

Primary [kms '] 

Secondary [km s * ] 

50080.5402 

-21.76 + 0.46 

30.06 + 0.72 

50081.5561 

-38.87 + 0.41 

50.28 + 0.64 

50090.6749 

-41.26 + 0.44 

51.80 + 0.68 

50094.6175 

56.08 + 0.66 

-54.20+ 1.04 

50098.5872 

-23.76 + 0.46 

32.20 + 0.72 

50108.5424 

-39.93 + 0.41 

49.52 + 0.64 

50114.6418 

41.45 + 0.43 

-38.19 + 0.67 

50115.4714 

7.68 + 0.43 

-1.97 + 0.67 

50118.4844 

-46.86 + 0.39 

56.77 + 0.62 

50120.5183 

-9.71+0.41 

18.45 + 0.65 

50121.5423 

52.49 + 0.42 

-50.82 + 0.66 

50126.6304 

-41.43 + 0.63 

50.12 + 0.99 

50127.5361 

-46.25 + 0.49 

58.01+0.77 

50140.5293 

83.47 + 0.41 

-81.43 + 0.65 

50143.5077 

-22.22 + 0.44 

32.63 + 0.68 

50144.5064 

-40.16 + 0.40 

51.78 + 0.63 

50146.5016 

-38.69 + 0.48 

49.69 + 0.75 

50147.5844 

-4.42 + 0.43 

12.59 + 0.67 

50152.4997 

-22.40 + 0.45 

33.31+0.71 

50153.5028 

-40.55 + 0.45 

51.12 + 0.70 

50154.5311 

-46.08 + 0.43 

58.03 + 0.67 

50155.4860 

-38.94 + 0.43 

48.45 + 0.68 

50156.4965 

-8.02 + 0.40 

16.87 + 0.63 

50170.5041 

-22.81+0.44 

33.28 + 0.68 

50173.5107 

-38.12 + 0.38 

48.41+0.60 

50176.5112 

83.02 + 0.44 

-82.74 + 0.69 

50179.5248 

-24.62 + 0.52 

34.14 + 0.82 

50336.8741 

25.92 + 0.42 

-22.09 + 0.66 

50346.8207 

81.42 + 0.42 

-80.38 + 0.66 

50348.7661 

23.78 + 0.45 

-17.67 + 0.70 

50350.8059 

-34.31+0.68 

44.49+ 1.07 

50352.8106 

-45.10 + 0.41 

54.89 + 0.65 

50356.8289 

64.68 + 0.42 

-62.93 + 0.66 

50360.7872 

-45.06 + 0.44 

55.20 + 0.68 

50363.8342 

26.13 + 0.43 

-21.03 + 0.68 

50443.6832 

-26.87 + 0.46 

36.25 + 0.73 

50449.5930 

-32.42 + 0.52 

41.67 + 0.82 

50456.5688 

26.22 + 0.49 

-20.28 + 0.77 

50460.5754 

-45.66 + 0.50 

56.81+0.78 

50462.4976 

13.02 + 0.58 

-5.35+0.91 

50464.6145 

67.50 + 0.50 

-66.16 + 0.79 

50466.6979 

-13.93 + 0.49 

21.60 + 0.77 

50472.5453 

77.67 + 0.50 

-76.34 + 0.78 

50474.6912 

20.80 + 0.50 

-15.25 + 0.79 

50477.5770 

-44.25 + 0.44 

55.52 + 0.69 

50478.5223 

-45.19 + 0.49 

56.62 + 0.77 

50481.6877 

81.27 + 0.59 

-82.90 + 0.93 

50493.5115 

-8.87 + 0.41 

15.61+0.65 

50495.4914 

-43.36 + 0.39 

55.08 + 0.62 

50503.5415 

-32.98 + 0.50 

43.68 + 0.78 

50505.5246 

-45.01 + 0.54 

55.96 + 0.84 

50516.5760 

21.87 + 0.56 

-17.99 + 0.87 

50523.5045 

-44.86 + 0.57 

56.30 + 0.89 

50527.5177 

68.51+0.57 

-66.85 + 0.89 

50531.5400 

-45.01 + 0.54 

54.87 + 0.84 

50535.4986 

77.96 + 0.55 

-76.90 + 0.86 

50536.5040 

68.85 + 0.54 

-67.56 + 0.84 

50538.5021 

-11.08 + 0.57 

17.92 + 0.89 

50540.5290 

-45.28 + 0.58 

55.64 + 0.91 
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3.2. Radial velocity measurements 

3.2.1. CfA spectroscopy 

All our CfA spectra appear double-lined. Radial velocities were 
obtaine d using the two-dimensi onal cross-correlation technique 
TODCOR (IZucker & Mazehlll994t) . with templates chosen from a 
large library of calc ulated spectra based on model atmosphere s 
by R. L. Kurucz Isee lNordstrom et al.ll994tlLatham et al.l2002h . 
The four main parameters of the templates are the effective 
temperature Tgff, projected equatorial rotational velocity vsini, 
metallicity [m/H], and surface gravity log^. The parameters that 
have the largest effect on the measured radial velocities are Tgff 
and vsini. Consequently, we held log^ fixed at 4.0 for the hot¬ 
ter and more massive star (hereafter star 1, or primary star) and 
4.5 for the cooler one (star 2, secondary star), which are near 
the final values reported below in Sect. 13.41 and we assumed so¬ 
lar metallicity. The optimum Tgff and v sin i values were deter¬ 
mined by running grids of cross-correlations, seeking the max¬ 
imum of the correlation coefficient averaged over a ll exposures 
and w eighted by the strength of each spectrum (see lTorres et alJ 
l2002h . The projected equatorial rotational velocities we obtained 
are u sin / = 9.0+2.0kms“' for star 1 and t; sin / = 4.4+2.0km s“' 
for star 2. The ysin/ value for star 2 is much less than the in¬ 
strumental broadening and so is very sensitive to the value of 
the macroturbulence assumed in the model grid. The effective 
temperatures derived are Tgffj = 5860 ± 100 K for star 1 and 
Tgff ,2 = 5780 ± lOOK for star 2. The uncertainty in these values 
has little effect on the measured radial velocities. 

As in previous studies using similar spectroscopic material, 
we made an assessment of potential systematic errors in our ra¬ 
dial velocities that may result from residual line blending as well 
as lines shifting in and out of our narrow spectr al window as a 
function of orbital phase (see lLatham et alJl996h . We did this by 
pe rforming numerical simulations analogous to those described 
by iTorres et"^ (Il997h . and we applied corrections to the raw 
velocities based on these simulations to mitigate the effect. The 
corrections were typically less than 0.5 km s ' for both stars. 

Finally, the stability of the zero-point of our velocity system 
was monitored by taking nightly exposures of the dusk and dawn 
sky, and small run-to-run corrections (typically under lkms~' ) 
were applied to the velocities as described by iLathaml (Il992h . 
The adopted heliocentric velocities including all corrections are 
listed in Tabled 

3.2.2. INT spectroscopy 

We also used todcor to measure the radial velocities of both 
stars from our INT spectra, but for these measurements we used 
spectra of HD 216435 (GOV) and HD 115617 ( G5V) obtained 
from a library of high-resolution stellar spectra (iBagnulo et alJ 
l2003h as templates for the primary and secondary stars, respec¬ 
tively. Both the INT spectra and the template spectra were inter¬ 
polated onto a uniform logarithmic wavelength grid equivalent 
to 7.23 km s ' per pixel using quadratic interpolation. The radial 
velocities were derived by fitting a minimum curvature surface 
to the peak of the two-dimensional cross correlation function and 
interpolating to the point of maximum correlation. The accuracy 
of the radial velocity measurements that can be obtained from 
our INT spectra is limited to about 1 kms^' by instrumental ef¬ 
fects (motion of the star in the slit, flexure, etc.) so we did not 
attempt to make any corrections for blending, etc. as we did for 
the CfA spectra. Radial velocities for spectra obtained near pri¬ 
mary or secondary eclipse were found to be unreliable and are 


Table 4. Radial velocity measurements from INT spectroscopy for 
VI094 Tau. The standard errors in the primary and secondary radial 
velocities estimated from the residuals for a least-squares fit of a Kep- 
lerian orhit are 1.4km s“' and 1.7 km s“', respectively. 


HJD-2400000 

Primary 

[km/s] 

Secondary 

[km/s] 

52561.5866 

-28.94 

32.28 

52561.5903 

-27.30 

33.68 

52561.5940 

-26.60 

34.94 

52561.5977 

-29.11 

32.46 

52561.6014 

-28.97 

32.56 

52561.7507 

-31.40 

35.83 

52561.7544 

-36.43 

31.43 

52561.7581 

-33.42 

34.26 

52561.7618 

-37.48 

30.59 

52561.7655 

-37.51 

30.45 

52561.7789 

-32.93 

35.66 

52561.7826 

-35.62 

32.86 

52562.5627 

-43.50 

49.21 

52562.5664 

-44.83 

47.26 

52562.5701 

-45.01 

47.22 

52562.6743 

-44.66 

50.30 

52562.6781 

-44.90 

50.23 

52562.6818 

-44.38 

51.03 

52562.7685 

-46.72 

50.16 

52562.7722 

-48.43 

48.51 

52562.7759 

-46.51 

50.33 

52563.6158 

-49.68 

53.96 

52563.6195 

-49.58 

54.07 

52563.6232 

-49.51 

53.58 

52563.7631 

-48.04 

54.07 

52563.7668 

-48.84 

52.88 

52563.7705 

-48.88 

52.98 

52566.7185 

67.66 

-71.19 

52566.7222 

66.88 

-72.05 

52566.7259 

66.96 

-72.84 

52566.7297 

67.19 

-72.72 

52566.7334 

67.76 

-72.67 

52566.7611 

69.18 

-74.22 

52566.7648 

69.41 

-74.14 

52566.7686 

70.05 

-73.39 

52566.7723 

70.28 

-74.14 

52566.7760 

70.54 

-74.16 

52566.7806 

70.42 

-74.19 

52566.7843 

69.21 

-76.27 

52566.7867 

70.69 

-77.47 

52568.6252 

33.26 

-35.50 

52568.6289 

33.09 

-34.80 

52568.6325 

32.77 

-34.66 

52568.7407 

28.53 

-30.12 

52568.7443 

27.30 

-31.13 

52568.7480 

26.81 

-31.41 

52570.7502 

-32.28 

35.55 

52570.7539 

-33.54 

34.33 


not reported here. The radial velocities derived from the 48 re¬ 
maining spectra are given in Table 0] 
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3.3. Spectroscopic orbit and apsida! motion 

We first performed least-squares fits of Keplerian orbits to the 
three sets of radial velocity data available to us, those from the 
CfA spectra, those fr om the INT/IDS s pectra , and the radial ve¬ 
locities published by iGriffin & BoffinI (i2003h . The results from 
the three data sets were found to be consistent with each other 
except for the differences in the zero-point of the radial veloc¬ 
ity scale between the different data sets. The offset between the 
INT radial velocities and the other radial velocity data is small 
compared to the instrumental resolution (14km s ') and is likely 
to be due to the uncertainty in the zero-point of the wavelength 
calibration for these spectra. The zero-point of the CfA radial ve¬ 
locity scale is within 0.14 km s ' of the absolute reference frame 
set by minor planets, which we have observed regularly with the 
same instrument for 25 years. 

There is clear evidence for apsidal motion in VI094 Tau 
from the measured times of mid-eclipse. There is also informa¬ 
tion about the rate of apsidal motion in the radial velocity data. 
In order to obtain the best possible estimate of the apsidal period, 
U, and to ensure that there is no systematic error in parameters 
of the spectroscopic orbit due to the variation in the longitude of 
periastron, w, we used a program called omdot to perform a si¬ 
multaneous least-squares fit to all three radial velocity data sets 
and all the measured times of mid-eclipse for the parameters of 
a Keplerian orbit in which the longitude of periastron changes 
at a constant rate at - InjU. The times of mid-eclipse are com¬ 
puted by calculating the times when the projected separation of 
the stars is at a minimum. The optim al solution is obta ined using 
the Levenburg-Marquardt algorithm (iPress et al.ll992h . We were 
careful to assign accurate standard errors to all the data so that 
the relative weighting of the different data sets is correct and the 
standard error estimates for the free parameters are accurate. The 
free parameters in the fit were: a reference time of mid-eclipse. 
To; the anomalistic period, Panom; the orbital eccentricty, e; the 
longitude of periastron at time To, wo; w; the semi-amplitudes of 
the spectroscopic orbits, K\ and K 2 ', the radial velocity of the bi¬ 
nary centre-of-mass for the CfA data, ycfA; two offsets between 
the different radial velocity scales, ycB -ycfA and yiNx-ycfA- The 
orbital inclination has a negligible effect on the derived times of 
minimum and the spectroscopic orbit so we fix this quantity at 
the value i = 88.25°. The sidereal period (mean time between 
eclipses) is PsW = (l/Panom + . The results from our least- 

squares fit are given in Table |5] 

The fit to the observed times of mid-eclipse is shown in 
Fig. 13 The variation in w over the time covered by radial ve¬ 
locity measurements turns out to be small (< 0 . 2 °) so we show 
the fit to the radial velocities for a fixed value of to and the best-fit 
values of Ki and K 2 in Fig. [2 

3.3.1. V-band and y-band luminosity ratios 

There can be some degeneracy between the ratio of the radii 
and the luminosity ratio determined from the light curve for par¬ 
tially eclipsing binaries, so it is useful to include the luminos¬ 
ity ratio from the spectroscopy as an additional constraint in the 
least-squares fit. The analysis of the CfA spectra using todcor 
yields a luminosity ratio in the 45A region around 5190A of 
^5190 = 0.534 + 0.020. To convert this value into a luminosity 
ratio in the V-ban d we used the synthe tic stellar spectra from the 
BaSel 3.1 library (IWestera et al.l2003) . We compare d the flux in¬ 
tegrated over the V-band response function (lBesselllll990l) to the 
flux integrated over the 45 A region around 5190A for a range of 
effective temperature and metallicity ([Fe/H]) similar to those 



-600 -400 -200 0 200 

Cycle number 


Fig. 2. Residuals from a linear ephemeris derived from the values of 
To and given in Table [5] (points) together with the fit to these val¬ 
ues from OMDOT for an eccentric orbit with apsidal motion (solid lines). 
Secondary eclipse is assumed to occur at phase 0.652 for this plot. Pub¬ 
lished times of mid-eclipse are plotted using diamond symbols, new 
times of mid-eclipse are plotted with dots (V-band) or open squares 
(SAT). 
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Phase from primary minimum 

Fig. 3. Radial velocities for the primary (filled symbols) and secondary 
(open symbols) components of VI094 Tau relative to their barycentre. 
Triangles(black) are Cf A data, circles (red) are INT data and squares 
(blue) are the data from IGriffin & BoffinI ( 120031) . Spectroscopic orbits 
are shown (solid lines) for parameters given in Table[5]but with to fixed 
at the value 333.61°. A colour version of this figure is available in the 
on-line version of this article. 


derived below for the two stars in VI094 Tau. By taking the 
ratio of these integrated fluxes and then the ratio of that ratio 
for the two stars we can derive a factor to convert ^5190 to a lu¬ 
minosity ratio in the V-band, fy. Applying this correction we 
obtain iy = 0.557 + 0.025, where the standard error estimate 


Article number, page 6 ofM 
































P. F. L. Maxted et al.: V1094 Tau 


Table 5. Parameters of the orbit for VI094 Tau derived using omdot. 
The number of eclipse times is A^tmin, the number of radial velocity mea¬ 
surements for each star are and Vrv ,2 and the standard deviation of 
the residuals for each of these data sets are (Ttmin, CTrv.i and cr„ 2 , respec¬ 
tively. The contribution of the eclipse times and radial velocity mea¬ 
surements for each star to the total value of for the least-squares fit 
aieXf^i^,Xni and^J^ 2 > respectively. Other symbols are defined in the 
text. 


Parameter 

Value 

Notes 

To 

2454537.54129 + 0.00019 


P anom [d] 

8.9885627 + 0.0000036 


e 

0.26767 + 0.00037 


coo n 

333.58 + 0.17 


d> [radians/period] 

0.0000107 + 0.0000027 


in 

88.25 

Eixed. 

Ki [kms-*] 

65.38 + 0.07 


fi:2[kms-'] 

70.83 + 0.12 


TCfA 

3.45 + 0.05 


Tgb - TCfA 

1.15 + 0.09 


Tint - TcfA 

-3.53 + 0.16 


Psid 

8.988547 


U[y] 

14500 + 3700 


Crtmin [d] 

0.00125 



41 


^ tmin 

39.2 


CTrv 1 [kms^'] 

0.93 


A^rv.l 

143 


U?v,l 

151.9 


crrv,2 [kms^'] 

1.17 


Nr.,2 

143 


U?v,2 

146.3 



337.4 



includes the additional uncertainty in this estimate due to the er¬ 
rors in the effective temperatures of the two stars and in [Fe/H]. 
A similar calculation for the Stromgren y-band gives the result 
{y = 0.533 + 0.025. If we were to extrapolate the luminosity ra¬ 
tio at 5190A to other wavelengths the results will be sensitive 
to systematic errors in the synthetic stellar spectra, so we have 
only used the luminosity ratio in the V-band and y-band using 
this method. 

3.3.2. Effective temperature and metallicity 

The two stars in VI094 Tau are quite similar and so a useful es¬ 
timate of the metallicity can be obtained by analysing the Strom- 
gremn photometry of this binary as though it wer e a single star . 
The reddening was estimated using the method of lOlse ^ (Il988h 
applied to the mean photometric colours out of eclipse from Ta¬ 
ble [T] and was found to be E(b - y) - 0.019 ± 0.007 magni¬ 
tudes. The calibration of iHolmberg et al.l d2007h applied to the 
de-reddened Stromgren photometry of VI094 Tau then yields 
the estimate [Fe/H] = -0.09 + 0.11, where the error includes both 
the random err or in the photometry an d the systematic error in 
the calibration. [Holmberg et alJ (l2007h also provide an effective 
temperature calibration which we can apply to our de-reddened 
Stromgremn photometry to obtain the value 5680 + 70 K. This 
value will be a weighted average of the individual effective tem¬ 
peratures for the two stars in VI094 Tau. The two stars have 


similar effective temperatures so the luminosity ratio is almost 
constant at optical wavelenghs, and so the relative weights be¬ 
tween the two stars will be close to the luminosity ratio derived 
above for the y-band. 

The surface brightness ratio derived from the least-squares 
fits to the light curves, J, provides a useful constraint on the 
difference between the effective temperatures of the two stars, 
ATfiff = Tsff.i - Tefr, 2 . We used the synthetic stellar spectra from 
iKurucj (Il993h to establish calibrations between surface bright¬ 
ness against Tes for the V-band, y-band and fe-band assuming 
either log o' = 4.0 (secondary) or \ogg - 4.5 (primary). As ex¬ 
pected, the value of log g has little effect on the predicted sur¬ 
face brightness. We then interpolated between these calibrations 
to find the value of Teff ,2 that gives the observed surface bright¬ 
ness ratio assuming Tesq = 5860 + 100 K. We investigated how 
these calibrations are affected by the assumed metal l icity u sing 
the synthetic stellar spectral library bv lWestera et al.l (l2002h and 
found that the estimate of ATefr changes by less than lOK if 
the assumed value of [Fe/H] is changed by +0.1 dex. There is 
very good agreement between the ATeff values derived using the 
two different spectral libraries. The values of ATefr derived by 
this method are 125 + 14 for the NFO V-band, 159 + 15 for the 
URSA V-band, 170+14 for the y-band, and 152 + 14 for the SAT 
b-band. The weighted mean of these estimates for the tempera¬ 
ture difference is ATgfr = 150 + lOK, where we have quoted the 
standard error in the weighted mean based on the scatter between 
the four input values. This is only slightly larger (by about 30%) 
than the standard error in the weighted mean calculated from the 
standard errors of these values. This shows that the level of sys¬ 
tematic error in our J values is low, i.e., the quoted errors on J 
are close to the true error in these values. The Stromgren u and v 
bands are strongly affected by details of the stellar models such 
as line blanketing and convection, so we have not attempted the 
same calculation at these wavelengths. 

The three constraints on the effective temperatures of the 
stars are shown in Fig. |4] together with our adopted values 
Teff.i = 5850 K and Tefr,2 = 5700 K. We have assumed in this 
plot that the value of Tefr derived from the Stromgren photome¬ 
try is a weighted average of the two individual Teff values with a 
ratio of weights equal to the luminosity ratio in the y-band. Given 
the level of agreement between the three constraints and the es¬ 
timated precision of the CfA effective temperature estimates, we 
have adopted estimated standard errors on both these values of 
+ 100 K. 

3.3.3. Light curve analysis 

We used jktebop version 25 (ISouthworthI 1201 3l and refer¬ 
ences therein) to analyse our six independent light curves of 
VI094 Tau (uyby, VNFn and Vtirsa) using the ebop light curve 
model (lEtze 1 1 198 It jPopper & Etzell 119811) . We only included 
data from nights covering the primary or secondary eclipse and 
we modified jktebop to include offsets in the zero-point of the 
magnitude scale between different nights as free parameters in 
the least-squares fit. This increases the number of free parame¬ 
ters in the least-squares fit, but it enables us to check that none of 
the parameters of interest is strongly affected by night-to-night 
offsets in the photometry. Eor light curves such as those pre¬ 
sented here where data from different nights overlap in phase 
we find that the offsets are well constrained and uncorrelated 
with other parameters so including them in the least-squares fits 
does not cause any problems. The offsets typically have values 
of a few milli-magnitudes with standard errors of about 1 milli- 
magnitude, though a few nights have offsets of about 10 milli- 
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their standard errors are given in Table |6] where we also pro¬ 
vide the weighted mean values and their standard errors for the 
wavelength-independent parameters. The standard errors quoted 
are the standard error in the weighted mean based on the scat¬ 
ter between the six input values. This is only slightly larger (by 
about 30%) than the standard error in the weighted mean cal¬ 
culated from the standard errors on the six input values. This 
shows that the level of systematic error in these free parame¬ 
ters is low, i.e., the quoted errors on the model parameters are 
close to the true error on these values. In particular, we checked 
that there is no trend in the estimated inclination versus wave¬ 
length that might suggest contamination of the light curves by 
third light with a different spectrum to the stars in VI094 Tau. 
The best fit to the V-band photometry from URSA around the 
eclipses is shown in Fig. |5] together with the residuals from the 
least-squares fits for all six light curves. 


Fig. 4. Adopted effective temperatures and standard errors (points with 
error bars) assuming [Fe/H] = -0.09. Constraints (l-cr) on the effective 
temperatures are indicated as follows: solid lines - CfA spectroscopy, 
dotted lines - Strbmgren photometry, dashed lines - V-band luminosity 
ratio. 


magnitudes. Other free parameters in the least-squares fit were: 
a normalisation constant, the surface brightness ratio J - S 2 IS 
where S 1 is the surface brightness of star 1 at the centre of the 
stellar disc and similarly for S 2 ; the sum of the radii relative to 
semi-major axis, {Ri + R 2 )la\ the ratio of the radii, k — R 2 IR 1 ', 
the orbital inclination, i; the phase of primary eclipse, A. We 
fixed the time of primary eclipse, the orbital period, e, and to by 
using the results from omdot described above to calculate these 
value of the quantities at the mid-point of the observed data. Sep¬ 
arate light curve solutions that included e cos(w) and e sin(fa;) 
as free parameters showed that the values derived are consistent 
with those derived using omdot. The variation in these quanti¬ 
ties during the span of the observations due to apsidal motion 
has a negligible effect on the results. For the V-band and y-band 
light curves the results presented here include the luminosity ra¬ 
tio calculated in section 13.3.11 as an additional constraint in the 
least-squares fit. We also tried least-squares fits without includ¬ 
ing the luminosity ratio as a constraint and found that this has a 
negligible effect on the parameters derived. 


Gravi ty darkening coeffic i ents, 0 , were taken from the tab¬ 
ulation of Iciaret & BloernenI (1201 lb , although gravity darken¬ 
ing has a negligible impact on the light curves for these nearly- 
spherical stars. We tried a quadratic limb darkening law for the 
fit to the b-bwd light curve and found that there was no im¬ 
provement in the quality of the fit compared to a linear limb 
darkening law and that the parameters of interest are negligi¬ 
bly affected, so we adopt linear limb darkening laws for all our 
light curves. The linear limb darkening coefficients for the two 
stars used in our analysis {xi and X 2 ) are given in Table|6] These 
are representative of the values obtained from several tabulations 
of limb-darkening coefficients for the approp riate effective tem¬ 
perature, gravity and metall i city of t he sta rs (vanjtonme 199^ 
Diaz-Cordoves et aH 119951 : IClarej I 2 OOOI: IClaret & Hauschii^ 


2003n . We estimate that the typical uncertainties on these val¬ 


ues are about 0.04. The uncertainty on the limb-darkening co¬ 
efficients has a small effect on the best-fit values of Ri/a and 
R 2 /a derived from each light curve (» 0.5%). We assume that 
these errors are random so that the effect on the weighted mean 
values of Ri/a and R 2 /a derived from the six light curves 
is negligible. The optimum values of the free parameters and 


3.4. Masses and radii 

We have used jktabsdiiv0 to combine the parameters of the spec¬ 
troscopic orbit from Table |5] and the weighted mean values of 
Ri/a, R 2 /a and i from Tabled to derive the masses and radii 
of the stars in VI094 Tau with their estimated standard errors 
given in Table]?] Table|7]also provides estimates of the luminos¬ 
ity and absolute V-band magnitude of the stars based on the es¬ 
timates of the stars’ effective temperatures in section l3'.3.2l The 
comparison with stellar models discussed below suggests that 
VI094 Tau may be slightly more metal-rich than assumed in 
section 13.3.21 The effective temperatures of the stars estimated 
from the CfA spectroscopy also increases by about 70K if the 
assumed values of [Fe/H] is increased by 0.1 dex. For that rea¬ 
son, we also provide in Table [T] estimates of the effective tem¬ 
perature, luminosity and absolute V-band magnitude of the stars 
for an assumed metallicity [Fe/H]=-i-0.14. In both cases, we also 
provide an estimate of the distance to VI094 Tau based on the 
surfa ce-brightness - effect ive temperature relation in the K-band 
from iKervella et af] (l2004ll . and the apparent K^-band magnitude 
of VI094 Tau from 2MASS transformed to the Johnson system 
(K = 7.468 + 0.021, LSkrutskie et al.ll2006HBesseiill2005h . 

4. Discussion 

Fig. ID shows the two stars in VI094 Tau in the mass- 
radius and mass-Te^ planes compared to various models from 
iBressan et al.l (120121) . The primary star is close to the end of its 
main-sequence lifetime and so models for this star are sensitive 
to the assumed age of the binary. For each value of the metal 
abundance shown we have adjusted the age of the models in 
order to find a good match to the radius of the more massive 
star. For a fixed metallicity we find that we can determine the 
age of the system with a precision of about 0.1 Gyr. If we as¬ 
sume, as Bressan et al. do, that the metal abundance of the Sun 
is Zq = 0.01774, and that the helium abundance of V1094 Tau is 
similar to the Sun, then our estimate [Fe/H]= -0.09 + 0.11 cor¬ 
responds to Z = 0.014 + 0.005. It is clear that for Z = 0.014 the 
effective temperatures of both stars are too low compared to the 
models. Increasing the assumed metal abunance to Z = 0.020 
provides a much better fit to the effective temperatures, particu¬ 
larly when we account for the increase in our estimates of Teg if 
the assumed metallicity is increased. A reasonable fit to all the 
observations can be obtained by assuming an intermediate value 
for the metal abundance, in which case the age of VI094 Tau is 

^ WWW.astro.keele.ac.uk/~jkt/codes/jktabsdim.html 
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Fig. 5. Upper panels: Detail of the jktebop model fit to the URSA V-band light curve of V1094 Tau around primary eclipse and secondary eclipse. 
Lower panels: Residuals from the jktebop model fit to our light curves of V1094 Tau, as labelled, in the region of the primary and secondary 
eclipse. For clarity, the residuals are offset in multiples of 0.1 magnitudes. Note that night-to-night offsets calculated as part of the least-squares 
fits have been applied to the observed magnitudes. 


Table 7. Absolute astrophysical parameters for both components of 
VI094 Tau determined wit h jktabsdim. Absol ute V magnitudes use 
bolometric corrections from iBessell et al.| lll998h . 



Primary 

Secondary 

Sidereal period [d] 

8.9885474: 

+ 0.000004 

Mass ratio 

0.9231 4 

: 0.0019 

Eccentricity 

0.26755 4 

: 0.00004 

Mass [Mq] 

1.0965 ± 0.0040 

1.0121 +0.0028 

Radius [Rq] 

1.4109 + 0.0058 

1.1063 + 0.0066 

log0 (cgs) 

4.179 + 0.004 

4.355 + 0.005 

Assuming [Fe/H] = 

-0.09 


Teff [K] 

5850+ 100 

5700 + 100 

log(L/L0) 

0.32 + 0.07 

0.05 + 0.03 

Mv 

4.00 + 0.09 

4.66 + 0.08 

Distance [pc] 

122 

+ 2 

Assuming [Fe/H] = 

+0.14 


Teff [K] 

5950 + 100 

5800 + 100 

log(L/L0) 

0.35 + 0.03 

0.10 + 0.03 

Mv 

3.91+0.09 

4.57 + 0.09 

Distance [pc] 

123 

+ 2 


estimated to be about 6 Gyr. Very s imilar results are found us ing 
other stellar mode ls grids such as IVandenBerg et al.l (l2006h or 
iDotter et al.l (l2008h . 

There is good agreement between the projected equatoria l 
velocity of the primary star quoted by iGriffin & BoffinI (l2003h 
(9 ± 1 km s ') and the value we have derived from our CfA spec¬ 
troscopy (9.0 ± 2.0kms '). The agreement is less good for the 
secondary star (7 + 1 km s“* versus 4.4+2.0km s“*) but the lower 
value derived from the CfA spectroscopy is very sensitive to the 


assumed macroturbulence parameter used in the stellar model 
grid because the rotational broadening in less than the resolution 
of the instrument (8.5kms“*). Nevertheless, both the primary 
and secondary stars appear to rotate slightly below the rate ex¬ 
pected for pseudo-synchronous rotation (11.42 ± 0.05 km s * and 
8.96 ± 0.05kms“', respectively). This suggests that V1094 Tau 
is close to the limit at which tidal interactions between the stars 
force them to rotate pseudo-synchronously with the orbital mo¬ 
tion. 


5. Conclusion 

We have measured the masses of the stars in VI094 Tau to bet¬ 
ter than 0.4% and the radii of these stars to better than 0.6%. 
The error estimates on these values are robust as they are based 
on the analysis of multiple high-quality independent data sets. 
This level of precision and accuracy is among the best currently 
available for any solar-type star apart from the Sun itself. The 
sample of 95 well -studied eclipsing binary stars compiled by 
iTorres et"^ (1201 Ol) contains only one star with both its mass and 
radius measured to this accuracy (HD 124784 B). 

We find that the current generation of stellar models are able 
to match the observed masses and radii of these stars within the 
constraints currently available on the effective temperatures of 
the stars and their metallicity. These models suggest that the 
age of VI094 Tau is about 6 Gyr. The rotation velocities of 
the stars suggest that VI094 Tau is close to the limit at which 
tidal interactions between the stars force them to rotate pseudo- 
synchronously with the orbit. 
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Table 6. Light curve parameters for V1094 Tau derived using jktebop. The number of observations used in the fit is N and the standard deviation 
of the residuals is cr. Other symbols are defined in the text. Values preceeded by “=” are fixed parameters or constraints in the least-squares fit. 



u 

V 

b 

L 

NFO 

URSA 

Adopted 

(Ri - 1 - R 2 )/a 

0.10804 

0.10830 

0.10870 

0.10870 

0.10737 

0.10772 

0.10805 

+ 

0.00068 

0.00043 

0.00047 

0.00035 

0.00028 

0.00058 

0.00026 

k = R2/R1 

0.751 

0.783 

0.799 

0.791 

0.788 

0.809 

0.788 

+ 

0.018 

0.014 

0.017 

0.010 

0.012 

0.018 

0.006 

R\la 

0.06169 

0.06075 

0.06044 

0.06069 

0.06005 

0.05954 

0.06050 

+ 

0.00051 

0.00036 

0.00043 

0.00028 

0.00032 

0.00048 

0.00024 

R 2 la 

0.04636 

0.04755 

0.04627 

0.04801 

0.04732 

0.04817 

0.04744 

+ 

0.00081 

0.00060 

0.00071 

0.00044 

0.00047 

0.00076 

0.00029 

e cos oj 

0.23969 

0.23956 

0.23957 

0.23962 

0.23985 

0.24000 


+ 

0.00012 

0.00007 

0.00007 

0.00006 

0.00004 

0.00010 


e sina> 

-0.11913 

-0.11915 

-0.11916 

-0.11924 

-0.11886 

-0.12079 


+ 

0.00003 

0.00004 

0.00004 

0.00006 

0.00060 

0.00095 


in 

88.31 

88.20 

88.20 

88.19 

88.22 

88.19 

88.210 

+ 

0.05 

0.03 

0.04 

0.03 

0.03 

0.04 

0.014 

J 

0.851 

0.850 

0.876 

0.899 

0.897 

0.873 


+ 

0.014 

0.008 

0.009 

0.007 

0.006 

0.009 


(= L2IU 

0.476 

0.517 

0.554 

0.559 

0.553 

0.567 


+ 

0.016 

0.014 

0.019 

0.011 

0.013 

0.020 


X\ 

= 0.821 

= 0.810 

= 0.753 

= 0.662 

= 0.662 

= 0.662 


X 2 

= 0.841 

= 0.822 

= 0.769 

= 0.679 

= 0.678 

= 0.678 


N 

492 

492 

492 

492 

3503 

5714 


cr [mtnag] 

8.5 

5.4 

5.3 

5.2 

6.5 

10.0 


Constraints 

e COS Cl) 

= 0.23970 

= 0.23970 

= 0.23970 

= 0.23970 

= 0.23976 

= 0.23964 


+ 

0.00031 

0.00031 

0.00031 

0.00031 

0.00045 

0.00045 


e sin Cl) 

= -0.11909 

= -0.11909 

= -0.11909 

= -0.11909 

= -0.11898 

= -0.11923 


+ 

{= L2/L1 

+ 

0.00039 

0.00039 

0.00039 

0.00039 

=0.533 

0.025 

0.00076 

=0.557 

0.025 

0.00076 

=0.557 

0.025 



fully in the data collection and analysis up to the time of his death, but bears no 
responsibility for the final text of this paper. 
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Fig. 6. Mass, radius and effective temperature of the two components of V1094 Tau (errror bars) compared to models from iBressan et all ( 120121) 
for various values of the initial metal abundance, Z. The assumed ages are 5.75 Gyr, 6.2 Gyr and 6.6 Gyr for Z = 0.014,0.017,0.020, respectively. 
The arrow in the lower-left panel indicates the effect of changing the assumed metallicity from [Fe/H]= -0.09 to [Fe/H]= -1-0.14. Small diamonds 
show the model grid points that have been interpolated to produce the isochrones in these plots. 
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